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Abstract 

Li showed that the R iemann Hypothesis is e quivalent to the nonnegativity of a certain 



sequence of numbers. iBombieri and Lagariad gave an arithmetic formula for the number 



sequence based on the Guinand-Weil explicit formula and showed that Li's criterion is equiv- 
alent to Weil's criterion for the Riemann Hypothesis. We provide a derivation of the explicit 
formula based on Laplace transforms and present an alternative expression for Li's criterion 
that invites a probabilistic interpretation. 
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1 Introduction 

The Laplace transform of a function f{t) on [0, oo) is dehned by 

m = C{f}{s) = / e-''f{t)dt 







If / converges for ?R.s = sq then it converges for all s with 5Rs > sq- 

Table [H sourced (directly or indirectly) from the more comprehensive table in [l|, p. 1020], 
provides Laplace transform pairs that are relevant to this paper. f*g denotes Laplace convolution 

{f*gm = f f{u)g{t-u)du 
Jo 

If f{t) is a real, non-negative function then F{t) = Jq f{u)du is called a distribution function 
with density f. F is a probability distribution and / the corresponding probability density if 
F{^) = 7(0) = 1. 

We take the view that the density / is the fundamental object of study because it allows the 
construction of a variety of (weighted) integrals involving / over prescribed intervals of [0, oo), 
the distribution F{t) being one such instance. This provides a simple, general and coherent 
framework for constructing explicit formulae of analytic number theory. 

The paper is largely expository in nature, using Laplace transforms to reproduce known 
results in analytic number theory with what might arguably be regarded as natural ease. The 
approach taken here is possibly best illustrated by example. Accordingly, we shall first study 
the Chebyshev counting function before proceeding to the general case. 
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Table 1: Laplace Transform pairs:- f{s) = e ^^f{t)dt 

2 The Chebyshev function: arithmetic form 

For 5Rs > 1, the Riemann zeta function is given by [gI, p. 6] 



n=l 



The rightmost form is the Euler product representation, where q runs over all the primes 
{2, 3, 5, 7, 11, . . .}. Its logarithmic derivative gives 



C'is) q-' f 1 , 



(2) 



g q 
With the aid of S4 in Table (TJ we recognise this as the Laplace transform of 

oo 

f{t) = Y,logqY,6{t-nlogq) t>0 (3) 

q n=l 

Hence / is a discrete density with an atom of strength logg at every multiple of logg. The 
associated distribution F is given by 

= ^logg^ / 6{u - nlog q)du = log q (4) 

The rightmost sum is a shorthand for a double sum over both q and n such that < e*. Setting 
X = e*, the Chebyshev counting function ip(x) is defined by 

ip{x) = F(logx) = log (7 (5) 

q"<x 

It is often written in terms of the von Mangoldt function A(n) as 



iP(x) = ^ A(n) A(n) 



n<x 



log q n = q^ 
otherwise 



(6) 
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In words, ip{x) is a distribution function over the integers with jumps of size log q at every power 
< X {k > 1). The density ([3]) can correspondingly be written as 

oo 

f{t) = Y,Mn)S{t- log n) t>0 (7) 

n=l 

For completeness, we note that ([2]) can also be written as 

C(s) ^ 

^ ' n=l 

The results presented here are well-known. The point to be noted is that we have arrived 
at V in a bottom-up fashion via the density /. We repeat the exercise below for the analytic 
representation of the zeta function. 

3 The Chebyshev function: analytic form 

The Riemann zeta and ^ functions are related by 0, p. 16]: 

2as) = s{s - l)vr-^/2r(./2)C(s) (9) 
Both ^ and F functions are expressible as infinite products 

oo ^ -I 



1-^ and -^^r(.) = n(l + 3" 

P n=l 



e 



where 7 is the Euler-Mascheroni constant and the product formula for ^ runs over all zeros {p} 
of 1^ (complex zeros of C) with p and 1 — p paired together. This formulation of the zeta function 
is valid for all s except for a simple pole at s = 1. 

Following the steps of the previous section, the logarithmic derivative of ^ gives 



— 1 ^-^ \s + 2n In I ^-^ s — p 1 

n=l ^ ^ p ' 



(10) 



p 

Laplace inversion leads to an alternative form for the density 



/(*) = ^* - E (^"'"* - ^ - E - ^(^ + i°g^)^(*) * ^ (11) 

n=l ^ ^ P 

Hence the corresponding distribution F is 

rt ^ —Int pt 11 



jgvrj 

n=l p p 



n=l p 



The second form follows from setting s = in (llOp . so that 



^'(0)- i + ^i_i(^ + iogvr) = -log2^ (13) 



C(0) ^ p 2 



3 



For a proof of the rightmost equahty, see [a, p. 67]. 
Finahy we may write ^l'{x) as 



%Ij{x) = F(log x) = X + ^ — ^ log 27r x > 1 (14) 

ZTl p 

n=l p 

This is known as von Mangoldt's expUcit formula for the Chebyshev function ip{x). The top- 
down approach to deducing the density from the distribution ^ or F is by differentiation 





oo 

= l_^^-2n-l 


-E 




n=l 


p 


or F'{t) 


oo 
n=l 


E=' 



^* t > 



While this may be valid for t > 0, the density remains unspecified at t = and therefore 
technically incomplete. But, as will emerge below, the atom at zero plays a significant role. 
To put the foregoing discussion in different words, we first note that, by G3 in Tabled 

Starting from F{s), we can 



1. Invert -F(sjj to obtain F{t) and then differentiate to obtain the density F'{t) valid for 
t > 0. 

2. Invert sF(s) to obtain the density directly. Here, we use G2 in Table [U 

£'^{sF{s)-F{+0)} = F'{t) =^ C-^{sF{s)} = F'{t) + F{+0)6{t) 

Hence the density we seek is F'{t) + F{+0)5{t), which is defined for all t > 0. This is 
equivalent to the density (fTT]l . with appropriate grouping of terms. 



Armed with a fully specified density, we now turn to the derivation of a more general explicit 
formula, subsuming von Mangoldt's explicit formula (|14p as a special case. 



4 The general explicit formula 

Given a function w{t) on [0, oo), Jq^ w{t)f{t)dt can be written in two ways using the arithmetic 
and analytic forms of the density f{t) given in ([3]) and pT]) respectively. We have already dealt 
in detail with the case w{t) = 1 above, where the integral yields the Chebyshev counting fuction 
in arithmetic and analytic form. 

In this section we shall primarily be interested in the limiting case T ^ oo. Subject to 
convergence, let {w{t)) denote the expectation of w{t) with respect to f{t) 

pT roo 

{w{t)) = lim / w{t)f{t)dt = / w{t)f{t)dt (15) 
"This is the Laplace transform equivalent of Edwards' explicit inversion [6], p. 50] to obtain il^i^) directly. 
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Consider w{t) = e ^^g{t) where g \s a. function on [0, oo) with Laplace transform g. The expec- 
tation with respect to the arithmetic density ^ gives 

{e-''g{t)) = hm /'°'" e-''g{t)f{t)dt (T = logx) 

/•oo 

= hm VA(n) / e'"^ g{t) 6 {t - log n)dt 

n<x ''^ 

•s—^ Kin) 

lim >^ — —g{\ogn) (16) 



X— >oo ^— ' 71 

n<a; 



The domain of vahdity in s wih depend on the choice of g{t)^ e.g. g{t) = e requires 5Rs > 1 — a. 
The expectation with respect to the analytic density pip gives 



f oo 

(e-^'git)) = I e-''g{t)f{t)dt 
Jo 

= e-^'g{t) l^e* - (e'^"* - i- 5{t)^ - ^ e^* - ^(7 + log 7r)5(t)^ 

°° / 1 \ 1 

= 5(5 - 1) - ^ U(s + 2n) - — 5(0) \ -J2g{s-p)--{^ + log 7r)5(0) (17) 

n=l ^ ' p 

This is the general explicit formula. The domain of validity in s will again depend on the choice 
of g(t\ although analytic continuation allows convergence in (jl7p for regions of s where (jl6p 
does not converge. We shall formally explore convergence for particular 'test functions' in the 
sections below, including the question of when we may legitimately equate (fT6]) and (fT7|) . For 
ease of exposition, we shall proceed with the general case for now without repeatedly calling to 
question its convergence properties. We shall have a specific interest in the case s = 1, so that 
(fTTll becomes 



^ 5(1 -P)= m - E (5(1 + 2n) - ^ 5(0)) - \ 

p n=l ^ ^ 



2(7 + logvr)5(0)-(e-*5(t)> (18) 



We proceed to derive an associated formula, using a related function g equivalent to what 
is referred to in [3] as the involution of g. We first note that G5 in Table [His normally defined 
for a > 0. It is also valid for a < provided that f{t/a) and f{as) remain meaningfully defined. 
For a = — 1, define 

g{t) = -e'g{-t) =^ g{s) = 5(1 - s) 



Hence 

/•oo 

{e-^'g{t)) = - / e~''g{-t)f{t)dt 



5(2-s)-^ U(l-s-2n) + — 5(0)j - ^ 5(1 " ^ + p) + ^ (7 + log 

n=l ^ ' p 



vr)5(0) (19) 



Setting s = 1 gives 



Y^~g{p) = 5(1) - f; (5(-2n) + i- 5(0)] + ^(7 + log7r)5(0) + {gi-t)) 

p n=l ^ ^ 



(20) 



Since ^ g{p) = ^ 5(1 — p) by virtue of the pairing of p and 1 — p, ([18]) and ([20]) are equivalent. 
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5 Weil's criterion 



We now turn to the convolution {g*g){t) whose Laplace transform is, by G6 in Tabled! g{s)g{s) 
g{l — s)g{s). Hence, noting that (g * g){0) = 0, the expectation of g * g( is 



r-OO 

{e-'\g*gm)= I e-'\g*g){t)f{t)dt 







oo 



= g{2 - s)g{s - 1) - ^ ^(1 - s - '^n)g{s + 2n) - ^(1 - s + p)g{s - p) (21) 

ra=l p 

Setting s = 1 as before gives 

oo 

Y,9{P)9{'^ -P)= 5(1)5(0) - Y.g{-2n)g{l + 2n) - (e-*(g * g){t)) (22) 

p n=l 

Weil's criterion states that a necessary and sufficient condition for the Riemann hypothesis to 
hold true is 

^5(P)5(1-P) >0 (23) 
p 

for all smooth functions g{t) [^@. 

It is clearly not feasible to test the criterion for all conceivable smooth functions. It is 
desirable therefore to find a manageable subset of test functions that can be shown to suffice. 
Li's criterion, discussed below, provides such a subset. By way of motivation, consider functions 
g that satisfy 

(g * gm = git) + g{t) = git) - e'gi-t) 

(24) 

=^ 9is)gis) = gis) + gis) = gis) + gil - s) 



Then the left side of (jlSp and (|20p can be written as 

E 9ip) = E - ^) = ^ E 9iP) +9il-p) = lYl 9iP)9il - P) (25) 



p p 



and (j22|l reduces to the sum of (jlSl) and (1201) . Hence, if a subclass of smooth functions satisfying 
()24|) can be shown to suffice for Weil's criterion, then the criterion amounts to the positivity of 
(fTSi) or, equivalently, (i20|) for such a set of functions. 

We now turn to some simple test functions to illustrate the ideas discussed thus far and to 
build toward the test functions needed for Li's criterion. 

6 Polynomial test functions 

The simplest case git) = 1 reproduces {e~^^git)) = /(s) = — C'('5)/C('5) which we already know 
can be written in both arithmetic dHj) and analytic form ()10p for 5Rs > 1 



s-l~ ^\s + 2n 2n) ^ s- p ^ 

71=1 n=l ^ ^ p ' 



(26) 



''The criterion generalises to complex smooth functions g{t) but real g{t) suffices for our purposes. 
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The choice g{t) = , k > gives {-)''f^^\s) 



n=l 



{s - 



'=ME 



+ E 



(«-p) 



fc+i 



(27) 



This shows that if g{t) is a constant, a power of t or a polynomial in t, then the explict formula 
amounts to working with derivatives of f{s) or a linear combination thereof for the polynomial 
case. This may seem to suggest that we can dispense with the explicit formula. However, the 
inadmissibility of s = 1 in ()26p and ()27p - where we have grouped together on the left side the 
terms that are not bounded at s = 1 - illustrates the benefit of retaining the view that the 
density / is the fundamental object of study. 

At s = 1 we retreat to the limit form (jlSp for the left side so that (j27p . say, becomes 



lim 



: lim 

X— >oo 



: lim 



n<x ■^'^ 

^A(n),, f^°s\ 

> logn y^- / t 

n<x 



s=l 



'dt 



E 

n<x 



A(n) 



n 



(logn)^ 



(log X 



,fc+i 



k + l 



(28) 



We have recognised the expression as (proportional to) the number % that arises in the Laurent 
expansiorQ about s = 1 of C,'{s)/C,{s) Hence (P7|) at s = 1 is 



°° 1 1 

^ (1 + 2n)k+^ ~ ^ (1 -p) 



^ (1 + 2n)'=+i 

n=l ^ ' 

It is straightforward to verify that 

^/ 1 1 



k+l 



k>0 



(29) 



^ (1 + 2 



^ (1 + 2n)^ 

n=l ^ ^ n=l 

Hence (j29p may be written as 



1= (l-2-MC(A:)-l k>l 



E(-i 



% + (-)^(l-2-'=-i)C(A; + !)-(-) 



> 



(30) 



The A: = case is already known from p3p to be 

= l + ^(7-log47r) 

P 

Finally, we consider g{t) = —e^g{—t) so that, for g{t) = t^, k > 



(31) 



^We avoid power series in s, in keeping with the following sentiment that Edwardj p. 9] attributes to 
Riemann: "The view of analytic continuation in terms of chains of disks and power series convergent in each disk 
descends from Weiersrass and is quite antithetical to Riemann's basic philosophy that analytic functions should 
be dealt with globally, not locally in terms of power series." 
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Hence, for /c > 0, (fT9]) becomes 



A; > (32) 



The s = 1 case does not require a limiting approach this time since all terms are bounded. Let 



Then 



becomes 



k+l 



k>0 



2-(fc+i)^(A; + i) _ k>0 



(33) 



(34) 



The A: = case simply reproduces (fSTI) . For completeness, we note that fiQ = — logvr. 



7 Li's Criterion 

Li 0] proved that a necessary and sufficient condition for the Riemann hypothesis to hold is 



A„ = ^l-(^l-i^ >0 n>0 



with p and 1 — p paired together. iBombieri and LagariasI [3] showed that Li's criterion is an 



instance of Weil's criterion restricted to a set of test functions {s'n(i)} defined as follows. 
The Laguerre polynomial L„(t) of degree n > and its Laplace transform are given by 



*;=0 ^ ^ 



i-t)^ t > 



1 

s / s \ s 



The associated Laguerre polynomial Qnit) = L]^_i{t), gn{t) = — e*L^_;^(— t) (n > 0) and corre- 
sponding Laplace transforms are defined bjQ 



n-l 
k=0 



k=0 



gn(s)=5n(l-s)=Ll„_l(l-s) =1-(^1-^) =1-(1 

Then, as observed in [3], gnis)gn{^ ~ s) = gnis) +fi'n(l — •?) by the identity (1 — r)(l — r~^) 
(1 — r) + (1 — r~^). Hence, as discussed above, Li' criterion is an instance of Weil's criterion 

An = '^9n{p) = ^^5n(p) +5n(l -P) = ^ J]] 5n (/?)?« (1 " p) > 

p p p 



''Also see [5|, iS] for identification of tlie test functions as tlie associated Laguerre polynomials Ll^_i{t). 



8 



To derive an explicit form, we use the expressions obtained in the previous section 

A.^E-(-i)"^-E(:)E(-i)'-Ei-E(:)E(-i 

= n + ^ (7 - log4vr) - X; h {vk^i - (-)^(l - 2-'^)C(^)) - ^ (fj 

k=2 ^ ^ k=2 ^ ^ 

which, with the aid of the identity J2k=o {'k)(-f = l-n + YJl=2 (fc)(-)^' = 0> leads to 
An = 1 + 5 (7 - log47r) + r ((-)'^(l - 2-'=)C(A;) - r?,_i 



(35) 



The corresponding exercise for the expressions involving {/i^} leads to 

An = 1 + ^ (7 - log 47r) + r {2-\{k) + ^u^.i) (36) 

k=2 ^ ^ 



Bombieri and LagariasI [3|] proved (I35p . It is an arithmetic interpretation of A^, involving as it 
does the zeta function at integer argument and {rjk} given by the arithmetic formula ()28|l4 

An arithmetic interpretation for the (formally equivalent) form (I36p is less direct because 
the sequence \u h\ (I33D a rises from analytic continuation rather than from an explicit arith- 
metic construct. lApostoll derived closed form expressions for C,^^\Q), from which closed form 
expressions for can be inferrec]^- 

A study of the properties of {//fc} might shed light on whether > for n > does 
indeed hold. If f{t) were a probability density then /^^ would be the k}'^ moment and therefore 
necessarily positive. This does not, of course, hold in this case because /(s) is the Laplace 
transform of the density f{t) only for s > 1. Nonetheless, a probabilistic interpretation arguably 
seems no less compelling than an arithmetic interpretation. 

8 Discussion 

It may seem natural to ask whether it makes sense for the analytic version (jlip of the density 
f{t) to have an atom at zero when it is meant to be equivalent to the arithmetic version ([3]), 
which is zero at t = since A(l) = (1 is not prime). 

In the conditioning notation of probability theory, it is more informative to write the arith- 
metic density as, say, f{t\q) - to be understood as "/(*) given the prime numbers as prior 
knowledge". Hence f{t\q) is discrete with atoms at the primes log q, where log 1 = is excluded 
by prior construction. Similarly, the analytic density may be written as f{t\p): ^^f{t) given the 
complex zeros p of the zeta function as prior knowledge". Such prior knowledge is defined in s 
space and it forces neither discreteness of f{t\p) nor zero mass at t = 0. 

'^It is shown in [sj, 1^ that the sequence {rik } can be generated recursively with the aid of the Stiehjes constants 

T'^-— ii^ 2^^^ kT^j 

bypassing, thereby, explicit dependence on the von Mangoldt function A(n). 

^Inevitably, {C'^^O)} ^-^id hence {/ifc} can also be written in terms of the Stieltjes constants. 
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While f{t\q) and f{t\p) are constructed in accordance with their associated prior assump- 
tions, ah integrals derived therefrom must necessarily be in agreement. It is, after all, the 
integrals such as the Chebyshev counting function that are the objects of ultimate interest. 

A deeper probabilistic approach to the Riemann Hypothesis, drawing specifically from the 
theory of infinitely divisible distributions, will be explored in a sequel to this paper. The work 
discussed here will be of direct relevance because the Laplace transform plays a central role in 
the treatment of infinitely divisible distributions on [0, oo) as discussed in |7l. llC|. 



9 Conclusion 

We have constructed arithmetic and analytic forms for the density on [0, oo) from which sums 
and integrals of interest may be surmised. We have shown that the general explicit formula 
follows effortlessly from the specification of such a density. We have further given an alternative 
form for Li's criterion. 
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